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Abstract. The accuracy of the measurements of some astrophysical dynamical systems allows 
to constrain the existence of incredibly small gravitational perturbations. In particular, the internal 
Solar System dynamics (planets, Earth-Moon) opens up the possibility, for the first time, to prove 
the abundance, mass and size, of dark sub-structures at the Earth vicinity. We find that adopting the 
standard dark matter density, its local distribution can be composed by sub-solar mass halos with 
no currently measurable dynamical consequences, regardless of the mini-halo fraction. On the other 
hand, it is possible to exclude the presence of dark streams with linear mass densities higher than 
/l st > 10~ 10 M©/AU (about the Earth mass spread along the diameter of the SS up to the Kuiper belt). 
In addition, we review the dynamics of wide binaries inside the dwarf spheroidal galaxies in the MW. 
The dynamics of such kind of binaries seem to be compatible with the presence of a huge fraction of 
dark sub-structure, thus their existence is not a sharp discriminant of the dark matter hypothesis as 
been claimed before. However, there are regimes where the constraints from different astrophysical 
systems may reveal the sub-structure mass function cut-off scale. 

Keywords: dark matter theory, power spectrum, particle physics-cosmology connection, dwarf 
galaxies 
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1 Introduction 



The A-CDM cosmological model assumes the existence of two unknown energy components; a cos- 
mological constant, A, and the Cold Dark Matter, CDM. This model is the paradigm on cosmology 
due to the high precision at which the observations of the large scale Universe are described (Cosmic 
Microwave Background radiation [1], large surveys of galaxies [2], gravitational lensing surveys [3], 
etc.). Such observations seem to be equally satisfied independently of the DM nature; although it has 
to be cold (low velocities), and it can have a weak interaction, at most, with the standard model of 
particles. 1 On the other hand, observations at sub-galactic scales are sensitive to the particle physics 
of the DM. 

In the current paradigm structure grows from bottom up; the first objects formed at high redshift 
were some small sub-galactic units that collapsed from small perturbations in the primordial density 
field. These clumps then merge hierarchically to form the large bound objects we observe in the local 
universe today. The properties, mass and scale, of the small structure is, to some extent, determined 
by the properties of the specific DM particle physics. Once a primordial power spectrum is fixed, 
its linear evolution leads to a mass power spectrum with a low end cut-off, at some wavenumber 
&f s (mainly determined by the particle free streaming (fs)). Therefore, one can associate a minimum 
mass, m m i n , for the smallest gravitationally bound structure that could be formed in the early universe, 
to the mass contained in a sphere of radius rf s ~ n/kf s . For WIMP's like particles, the mass cut-off 
is set in the range of ~ 10~ 4 - 1(T 12 Mq [9, 10], while for Axions is set around 10~ 18 - 
10 Mq [11, 12], (the specific values depend on the decoupling temperature of the specific model). 

The structure formation process makes very difficult to confront the galactic and sub-galactic 
scales, which are non-linear, with observations. However, these scales are the ones that contains the 
more valuable information. To confirm the existence of mini-halos (sub-halos with a mass as small 
as predicted by the power spectrum cut-off), is crucial for the determination of the DM nature, not 
only because it could restrict the particle physics of DM, but also because in the case they do exist, 
they could affect the ongoing direct and indirect searches for DM. 

Numerical simulations have been very useful tools to study the non-linear evolution of structure 
in the larger scales. But, the survival and distribution of sub-halos (not only the mini-halos but 
in general), are difficult to study under this scheme, due to resolution limits and the difficulty in 
including the relevant baryonic process [13-15]. Also, semi-analytic formalisms have been used to 
study the assembly history of halos, alone or in combination with N-Body simulations, [16, 17]. 
These formalisms treats the sub-halos as isolated systems, not taking into account they were part 
of a bigger structure. Numerical simulations and semi-analytical formalisms seem to be in good 
agreement, but on scales much larger than the ones we are interested in this work. 
Also, some of the observational efforts, planed and/or under development, to constrain the mass 
power spectrum at galactic and sub-galactic scales are: fits to the fluctuations of gravitational lensing 
observables, like the time delays and the cusp-fold relation [18, 19], and the identification of ultra faint 
galaxies in current, and future, surveys [20, 21]. None of the above would actually probe the size of 
the mass cut-off. Perhaps the one suggested in [22] where the authors show there is a correlation 
between the mass cut-off and the direct detection rates, and/or fluxes, of high energy neutrinos from 
DM annihilation in the Sun; or the one in [23], where a correlation between the mass cut-off and the 
radius, and surface density, of the galactic halo is proposed, are the observational prospects, up today, 
that would probe the smallest scales of the power spectrum and their relation with the DM nature. 

'Particles like the sterile neutrinos can play a Warm DM role, which is still viable according to different astrophys- 
ical constraints [4, 5], The experimental evidence of neutrino oscillations admits its presence, though the cosmological 
constraints are not conclusive yet, see [6-8] for detailed analyses. 
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All of the above has motivated several studies about the evolution, specifically the survival, 
of the mini-halos under the dynamical interaction with stars, the galactic potential, or bigger sub- 
halos [24-26]. Some of these studies reveal that a huge fraction of the mini-halos survive the as- 
sembling process [27, 28], while others suggest that only a very small fraction of them remain as 
mini-halos. Nonetheless, the disruption process leads to stream like structures as remnants [29, 30]. 

In the present work, we go through a complementary methodology that allows to set limits to 
the mass, and distribution of mini-halos, or streams, in the Galactic halo. Previous studies have focus 
on the disruption of the mini-halos (or, the formation of the streams), by their interaction with stars. 
Instead, we will focus on what happens to the dynamical system bound to these stars, like planetary 
systems, binary companions, etc., due to the dynamical interaction with the sub-structure, — mini- 
halos, or already formed streams — . In general, we will study the dynamical stability of astrophysical 
systems by considering the rate of injection energy due to the tidal interaction with a family of DM 
mini-halos and/or streams (perturbers), whose size and abundance can be related to properties of 
some DM candidates. 

Our own planetary system is a good candidate to explore this approach, since the orbital ele- 
ments are known with outstanding accuracy. For instance, the Earth-Moon distance has been mea- 
sured with increasing precision, thanks to the APOLLO-LLR mission [31], such that has been used 
to test gravity theories for a long time. We will estimate the dynamical perturbations that inner Solar 
System dynamics, Earth-Moon and Sun-Neptune systems specifically, could suffer by the interaction 
with a background population of dark sub-structure; and compare with the current accuracy threshold 
in order to obtain constraints on the sub-structure properties. Other good systems to use within this 
approach are the wide binaries inside dwarf Milky Way (MW) satellites. The dynamical interaction 
between such kind of binaries and DM has been discussed in [32] and [33], and their survival has 
been proposed as a discriminant of the DM hypothesis. Here we review the dynamics of this systems 
under the dynamical interaction with mini-halo mass, and/or streams, and study the effect of multiple 
encounters. 

This work is structured as follows: first, we present the basic equations to compute the energy 
perturbation on an astrophysical system due to a single encounter with a mini-halo like perturber 
(section 2.1), and with a stream like perturber (section 2.2). Then we describe the treatment of 
multiple encounters by means of Monte-Carlo experiments in section 2.3. In section 3.1, and 3.2, 
we use the described method to study the dynamical perturbations to the inner Solar System, and to 
the open binaries, respectively. In section 4 we outline a possible connection between the linear mass 
density of streams, with the mass of the initial mini-halo. Finally, we discuss our results in section 5. 

2 Dynamic perturbations on astrophysical systems due to dark sub-structure 

The dynamics of some astrophysical systems is known with great accuracy, such as our own So- 
lar System, that any dynamical perturbation that could affect it could be constrained. Including the 
presence of dark sub-structure. These systems would experience an almost negligible effect if DM 
were smoothly distributed over the dark halo, however, if DM were mostly in form of mini-halos, or 
streams, these act as perturbers to the astrophysical system (the target). The effect of one single en- 
counter could be small, but the cumulative effect of all the halo sub-structure could play a significant 
role to the target dynamics. Below, we work out a formalism that allow us to quantify this effect for 
mini-halos and streams acting as perturbers. We will work with binary systems as targets, but the 
reader should keep in mind that all of what is presented here could be applied to any astrophysical 
system, just by describing it in the center of mass reference frame. 
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2.1 Single encounters with mini-halos 

If mini-halos has not suffered important disruption, they can be treated as spherical perturbers. For 
encounters that occur at distances much larger than the characteristic size of the target, and when 
the effective interaction time is a lot shorter than the internal dynamical time of the target,- say the 
orbital period of a planet-, the encounter is said to be in the distant and impulsive regime. Within 
this regime, the average energy perturbation, per unit mass, to an extended system due to one single 
encounter is: 

7G 2 m 2 . „ 

(AE)mh = - 2 f- a 2 U(pfr h ). (2.1) 
3vf )P 4 

where, the impact parameter p is the shortest distance from the target center to the mini-halo, is 
the mass of the mini-halo, vo is the relative encounter velocity, a is a characteristic size of the target 
(e.g. the mean particle separation), and r/, the characteristic radius of the perturber [34]. The U func- 
tion is a dimensionless correction factor that takes into account finite-size effects of the perturber (eq. 
(3) in [35]); it depends on the density profile of the perturber and takes values between and 1. As 
the impact parameter becomes larger than the characteristic size of the perturber, the U function tends 
to one and the point mass approximation becomes valid. Most of the encounters in our study falls 
within this approximation, for the rest we use the U function given in appendix A. 1 . The fractional 
energy perturbation to the target is the average energy perturbation normalized to the binding energy 
(|i?bl = GM c /(2a), with M c the mass of the central body in the target), given by: 

\E h \ 3M c vlp* a - 1 } 

If it happens that the distant approximation does not hold, then we implement the correction described 
in appendix A.2. 

2.2 Single encounters with Streams 

Streams are former DM mini-halos (sub-halos in general), that have been already disrupted by tidal 
forces. The interaction of streams with a bound system of finite size can be seen as the gravitational 
interaction of it with a structure of cylindrical symmetry. The simplest stream one can model is 1- 
dimensional (ID), i.e. a line of longitude L and linear mass density A st . Nevertheless more realistic 
models include: the finite cross section stream with constant density (CD); and the finite cross section 
with a central core (Core) or with a central cusp (Cusp). In equations (2.3) we present the equations 
we will use to compute the fractional energy perturbation to the target (normalized to its binding 
energy), for the different stream models. 
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(2.3a) 
(2.3b) 
(2.3c) 
(2.3d) 



The functions T, B, C and D are dimensionless and contains information about the geometry of the 
encounter, (defined by the angles 6, a, (3 and ip), and about the characteristic radius of the stream, 
Ro- They are the equivalent of the [/-function. All of the equations for the finite cross-section 
streams converge to that of the 1 -dimensional stream for large impact parameters. The reader is 
referred to appendix B. 1 for a derivation of equations (2.3) and for a wider explanation of the different 
parameters. 
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2.3 Multiple Encounters: Monte-Carlo simulations of the encounters 

Regardless of whether the interaction is with mini-halos or with streams, the fractional energy of the 
target system will change as a result of repeated encounters. Here we can distinguish 2 contrasting 
regimes: in the first, the "diffusive" regime, the binding energy slowly, but steadily changes as the 
effect of small individual encounter perturbations accumulates. In the second, the "catastrophic" 
regime, a few encounters are enough to unbind the target. Our general approach to study the effect 
of multiple encounters consists of Monte-Carlo experiments, which automatically accounts for the 
catastrophic and diffusive regimes. These are performed as follows: 

1. First, we construct individual histories of encounters. 

(a) We start bay drawing random values for the encounter parameters, (impact parameter, 
encounter velocity and orientation of the target, etc.), following their corresponding dis- 
tribution functions. 

(b) We use these parameters to assess the effect of one single encounter trough the fractional 
energy given by equation (2.2) for mini-halos, or equations (2.3) if the perturber is a 
stream. 

(c) Now we update the mean size of the target according to the fractional energy injection 
just obtained. The relation between fractional energy input, and the changes in the mean 
separation of the target is \AE/E\ = \Aa/a\ (from the definition, and the derivative, of the 
binding energy). 

(d) Then we increment the experiment timer by approximately the time the perturber takes to 
traverse a sphere of radius equal to the impact parameter t, with a encounter velocity vo: 

Afenc = 2p/v . 

(e) At last we repeat steps (a) to (d) until the time limit (e.g. the life time of the target), is 
reached, or until the fractional energy reaches the unity (i.e. the target gets unbound). 

2. The above will give us the total fractional perturbation energy to the target, due to one single 
history of multiple encounters between our target and a given realization of the perturbers 
distribution. The second part of our Monte-Carlo consist of creating an ensemble of such 
histories, (i.e. to repeat the above procedure a number of times AO, so we get to know the mean 
value of the total injection energy that a population of sub-structure produces on the dynamics 
of a specific target. 

Distribution of encounter parameters. 

Let us describe now the distribution functions we use to sample the different parameters. For encoun- 
ters with a mini-halos, the impact parameter will be sampled from a nearest-neighbor distribution, 
(random, Poisson distribution), modulated by the DM local density (pdmX the mass of the perturbers 
(m m h) and the fraction (f s ) of DM density that is in form of mini-halos: 

/ (p) dp - 4nv p 2 exp 3 , with v = s ^ dm . (2.4) 

"Imh 

On the other hand, for an encounter with a stream we sample the impact parameter from the proba- 
bility distribution function shown in Figure 13b. We can not use Eq.(2.4), now modulated now by the 
linear density of streams, because the final impact parameter distribution does depend on the sample 
volume used to generate it, and there is no way to set it a priori. Instead, we adopted a more basic 
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Figure 1: Resultant probability distribution of fractional energy imprinted to the Earth-Moon system 
due to multiple encounters with: (a) mini-halos of mass - 1O~ 6 M0/AU, and (b) streams of 
linear density A st = 10~ 6 Mq/A\J. The ensemble consists of about 10,000 histories of the target 
evolution, under the interaction with the corresponding perturber; for about 4 Gyr each history. 



procedure that does not depend on any sample volume, it is described in appendix B.2. 

For both types of encounters, with mini-halos or streams, we sample the relative encounter 
velocity from a Maxwell-Boltzmann distribution function. We are considering that the distribution 
of sub-structure is relaxed and well described by a classic distribution [30]. Finally, the encounter 
geometry parameters (6, a, etc.) are sampled from uniform probability distributions. 

3 Application to astrophysical systems 

If a population of DM mini-halos and/or streams is pervasive through the Galaxy, as some variants 
of DM constituents suggest, it is likely that their presence could be inferred indirectly from their 
cumulative dynamical effect on some astrophysical systems, which could be very sensitive to external 
perturbations. We now explore the consequences for two such systems. 

3.1 Solar System: Earth-Moon and Neptune orbits 

The Solar System dynamics is known with high accuracy. For instance, the uncertainty in the planet 
semi-major axis is of the order of 10~ 9 AU for Jupiter and 10~ 6 AU for Neptune [36] while for the 
Earth-Moon system, it is currently determined with a precision of Aa/a ~ 10~ 13 [31]. This precision 
has made the Solar System to be extensively used to test general relativity, fundamental physics and 
alternative theories of gravity (see [37, 38, and many others] for recent examples), and it is particularly 
attractive to constraint the amount of DM that may be present in the Solar System [39^12]. It is 
interesting as well to constrain the presence of dark sub-structure in the solar neighborhood. 

As the Solar System moves within the Galactic halo, the orbits of planets in this system will 
experience almost negligible perturbations if DM is smoothly distributed, contrary to the case if it is 
distributed in mini-halos or streams. In the later case, the planet orbits will undergo an energy pertur- 
bation due to multiple encounters with these structures, causing them to undergo secular evolution, 
and even to get unbound. We will compare the predicted dynamical effect, changes in the semi-major 
axis of planet orbits, principally, with the current accuracy threshold to determine upper limits on the 
DM sub-structure that could be present. 
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Figure 2: Constraints to the mini-halo mass due to encounters with the (a) Earth-Moon and (b) Nep- 
tune System. The median energy perturbation to the orbit, due to interactions with mini-halos as a 
function of their mass, compared to the experimental uncertainty in the semi-major axis of the corre- 
sponding orbit[31, 36] (dashed line). Three different values of the local DM density, pdmCMq/AU 3 ), 
are shown for comparison. The bars correspond to the interquartile range for the Montecarlo real- 
ization. The sub-panel shows the effect of include the correction for extended perturbers, and that of 
varying the fraction of sub-structure. 



To start with, we choose the Sun-Neptune and Earth-Moon relative orbits as targets. Both have 
very stringent constraints in the determination of their semi-major axis, and they are also suitable to 
work within the impulsive encounter approximation. Assuming that the encounter velocity is roughly 
200km/s, (the Sun goes around the galactic center at about 220km/s and the dark halo presumably 
has very little rotation), means that a perturber crosses Neptune's orbit in just 1.4 years, that is many 
times shorter than Neptune's orbital period of 164.8 years; we are definitively in the impulsive regime, 
a simmilar argument applies for the Earth-Moon System. When considering multiple encounters, 
we assume the encounter relative velocity magnitude to follow a Maxwell-Boltzmann distribution 
with a mean of 220km/s, as it is usual. Finally, we assume a standard local DM density value, 
p dm = 10" 18 M© AIT 3 (0.35 GeV/cm 3 ) [43-47] (unless otherwise is specified). 

In Figure (1) we show a couple of examples of the resultant distribution for the total fractional 
energy, AE/Eb, obtained in our Monte-Carlo experiments ensembles. Each of the runs in the en- 
semble corresponds to a history of encounters between the Earth-Moon system and a population of 
mini-halos of m mh = 10~ 6 Mo/AU (Figure ( la)), or streams of A st = 10 _6 Mo/AU (Figure (lb)). 
The evolution time of the system in each history is for about 4 Gyr, approximately the age of the 
Earth-Moon system, and the time that Neptune has been in its present orbit [48, 49]. The ensembles 
consists of about 10000 histories of the target evolution. For both systems the resultant distribution 
is clearly non-symmetric, then we will further report the results in terms of the median value and the 
interquartile range. 

First, lets take a look at the results for encounters with mini-halos. In Figures (2a), and (2b), 
we summarize them, for the Earth-Moon and Neptune-sun systems, respectively. The perturber mass 
is in the range 10 < m m h < 10 . We can see that none of both systems are sensitive enough to 
the presence of this kind of sub-structure, neither at the standard local DM density (blue line in the 
figures), nor at a greater values (red and green lines). Except for the most massive perturbers, more 
than 10 2 M©, and in the denser case, pd m > 10~ 18 M©. This would be the maximum effect that DM 
mini-halos can produce, since we have assumed the fraction of mass in mini-halos to be the unity 



-7- 



-4 








fsMU 



_ 16 L. — . — . — . — . — . — . — . — . — . — . — . — i — , — . — . — i — , — . — . — . — . — . — . — iJ Lt_^ — . — . — , — . — . — . — . — . — . — . — . — . — . — . — . — . — . — . — , — . — . — , — J 

-12 -10 -8 -6 -4 -2 -12 -10 -8 -6 -4 -2 

A st [M /AU] A st [Ms/AU] 

(a) Earth-Moon (b) Neptune 

Figure 3: Constraints to linear density of dark streams due to encounters with the Earth-Moon Sys- 
tem: We show the median energy perturbation to (a) the Earth-Moon's orbit and (b) Neptune's Orbit, 
due to interactions with dark mini-halos as a function of its mass and for three different values of 
the local DM density, pdm(M0/AU 3 ). The bars correspond to the first interquartile range for the 
Montecarlo realization, and the dashed line to the experimental uncertainty in the semi-major axis of 
the orbit. 



and treated them as point masses. To consider extended mini-halos, and smaller fractions of sub- 
structure, would result in lower values to the injection energy, (see appendix A.l, and the sub-panel 
in Figure (2b)). 

Now, lets turn our attention to the interaction of these systems with dark streams; the linear 
mass density values used here are in the range 10~ 12 < /t st (M©/AU) < 10°. Figures (3a), and (3b), 
resumes the results. Contrary to the point mass case, both systems could be significantly affected 
by the presence of dark streams. At the standard value of the local DM density (blue line in the 
figures), the Earth-Moon system excludes the possibility of having the DM distributed on streams 
with A st > 10~ 10 Mo/AU, while the Sun-Neptune does not impose any restriction. On the other hand, 
if we consider higher values for the local DM density (red and green lines in figures) the constraints 
becomes more stringent. Both systems will be compatible with the presence of streams with linear 
density A st < KT 1 'Mq/AU, in the denser case. We initially assumed all the DM density is in streams, 
f s = 1, once we relax this condition, (shaded blue area in the plots), there are several combinations 
of sub-structure fraction and stream linear densities that are allowed by the observational constraints. 
The allowed and not allowed combinations are shown in Figure(4): a density plot for the difference 
between the calculated fractional change to the semi-major axis of the orbit and the observational 
restriction to the same quantity, the solid line indicates where the two quantities are equal, separating 
the region of the parameter space that is excluded (yellow-red) from that allowed (green-blue). 

A more realistic situation would be one in which a fraction of the mass is in the form of streams 
while the rest is in mini-halos. Since the encounters with one or the other is independent. Then, in 
general the results will be set by the dominant specie. We can see from the plots that the dominant 
contribution would come from the streams in most of the cases. 

Evolution of the Astronomical Unit There is evidence of some unexplained anomalies that con- 
cerns the astrometric data, including: the Pioneer anomaly, the variations in the astronomical unit, the 
Earth flyby anomaly, and the increase in the eccentricity of the Moon's orbit (all of them discussed 
in [50]) . The first one has been already explained by reanalyzing the thermal radiation pressure 
forces inherent to the spacecraft [51], thus it has been discarded the possibility that the presence of 
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Figure 4: Fractional injection energy, for the Earth-Moon system, as a function of f s and A st . The 
solid line is where the perturbation equals the observational limit, assuming the standard DM local 
density. The right axis shows a possible relation with the initial sub-halo mass, see section 4. 

some dark matter is the responsible for this anomaly. The second one remains unexplained, (see 
however [52]). The increase in the astronomical unit is reported to be about 15 cm yr~ l [50]. We 
have applied our formalism of multiple encounters with dark sub-structure to the Sun-Earth system, 
to check if this anomaly could be explained by the dynamic perturbations due to the presence of sub- 
structure. Assuming the standard local dark matter density value, and regardless the sub-structure 
fraction, neither the perturbations due to mini-halos nor due to streams, for aboutthe lifetime of the 
Earth, would be energetic enough to explain this anomaly. If larger values of the local density were 
allowed, the anomaly could be explained, however it would cause the secular evolution of other sys- 
tems, like the Earth-Moon. 

3.2 Wide binaries in dSphs Galaxies 

It is known that binary stars systems with separations as large as ~ 0. 1 pc exist in the stellar halo of 
the MW [53, 54], but corroborating their existence in dwarf spheroidal galaxies is an unachieved goal. 
The survival of these wide binaries against the interaction with the DM has been discussed before. 
Those in the MW had been used to set strong limits to the MACHO's, (masive compact halos), 
mass[55]. And now, the existence of those in the dSph's had been proposed as a discriminant for the 
dark matter hypotesis [32, 33]. In particular in [32] it is argued that such kind of wide binaries should 
be disrupted due to interaction with DM sub-halos. But, only catastrophic encounters, within the 
tidal interaction approach, were considered. They also assume a sub-halo mass function like the one 
founded in Af-body simulations. It is important to recall that the formation process of wide binaries, 
and the determination of their lifetime is unclear [56, 57], as it is the assembling history of the dSph's. 
This uncertainties could play an important role, when estimating the dynamical perturbations that the 
binaries suffer, due to the presence of dark sub-structure. According to previous studies, the existence 
of wide open binaries in dSphs would not be compatible with the DM hypothesis. 

We study the disruption of wide binaries by means of their tidal interactions with two types of 
DM sub-structure: mini-halos and streams. Considering only the effect of sub-structure that is bound 
to the dwarf halo, and not that bounded to the MW. At the typical distance were the dSph's galaxies 
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Figure 5: The median energy perturbation to a binary star system due to interaction with mini-halos. 

(a) As a function of the mini-halo mass, m^M©) , and for three different evolution times, ? e v(Gyr). 

(b) As a function of the fraction, f s , of the mass density in mini-halos of different mass. 



are located, the density of the MW halo has drop enough to be taken into account. We will follow the 
cumulative effect of interactions with sub-structures by means of our Monte-Carlo experiments, this 
will automatically take into account the diffusive and the catastrophic regime. 

Such binaries could be present at a radius r ~ O.lkpc from the center of the dSph galaxy, 
where a typical value for the DM density is p^m ~ 2.2 * 10~ 17 M©/AU 3 , (assuming the galaxy is 
immersed in a NFW density profile with virial mass M v ; r = 10 9 M© and concentration c = 23.1, as 
in [32]). 2 Dwarf galaxies have typical dispersion velocities of cr « lOkm/s [58]. We then assume 
the relative velocity of encounters to follow a Maxwell-Boltzmann velocity distribution, as in the 
previous section, but with a dispersion of V2o~ ~ 14 km/s. Below, we present the results obtained by 
following the evolution of the binaries from their interaction with dark sub-structure: 

Interaction with mini-halos In Figure (5a) we report the internal energy change, of the binary sys- 
tem, as a function of the perturber mass in the range 1(T 4 < m m h < 10 8 M©. For three different 
values of the evolution time: the solid line corresponds to the typical value assumed in previous 
works (t ev ~ lOGyr), and the dashed, and dotted lines correspond to shorter times (5Gyr and 
2Gyr). The line at 5E/Eb = represents the limit at which the injection energy due to encoun- 
ters equals the binding energy of the binary, i.e. the moment at which the binary gets unbound. 
We can see that imposing a constrain to the lifetime of the binaries would impose a limit to the 
maximum mass of the sub-structure that could be present in the dSph's. However, this mass 
limit will also depend on the fraction of the mass that is in form of mini-halos. In Figure 5b we 
show the effect of varying the sub-structure fraction for selected values of the mini-halo mass. 
Note that in the case of having massive perturbers, > 10 2 M© (the mass limit between the dif- 
fusive and catastrophic regime for tidal encounters set in [32]) it would be possible to unbound 
the binary, only if all the DM were distributed in mini-halos of this mass, and the evolution 
time is as large as 10 Gyr. The probability of having a massive perturber at short distance to the 
binary decrease as the impact parameter has to be larger than for the less massive perturbers. 

2 This particular radius is smaller, in general, that the half-light radius of dwarf galaxies [58], and that the one assumed 
in [32]. We choose it this way because at smaller radius, there are more stars, but also, the DM density is higher, so the 
interactions could be more frequent. 




Figure 6: The median energy perturbation to a binary star system due to interaction with dark streams. 

(a) As a function of the stream linear mass density A st and for three different evolution times , ? ev (Gyr). 

(b) As a function of the fraction of the mass density in streams, for four different linear mass density 
values. 



Then it would be more probable to cause the disruption of the binary by the presence of small 
mini-halos, than having a catastrophic encounter. 

Interaction with streams For the interaction with streams we do the same as above, but, now as a 
function of the linear mass density A st . The results are shown in Figures (6a) and (6b). Again, 
the binaries could be dissolved for several combinations of A sU t ev and f s . Limits to the lifetime 
of the binaries could impose some limits to the fraction of the mass that can be composed of 
dark streams of a particular linear density. For instance, if we consider that binaries have been 
for about 2Gyr under disruption, then there can not be present streams with linear densities 
greater that 1(T 6 M0/AU, composing more that the 10% of the total dark matter density, oth- 
erwise the binaries would be disrupted. Similarly, if the binaries have been in the perturbation 
field for about lOGyr, then there could be present streams with A st < 1O~ 1O M0/AU, or even 
Am < 1O -8 M0/AU, depending on value of the / s parameter. It is interesting to note that the 
limit to the linear mass density of streams in the later case is quite similar to that imposed by 
the Earth- Moon system in section 3.1. 

A more general situation is that in which a fraction of the DM density mass is distributed along 
streams, and other in mini-halos. In this case the total energy perturbation will depend on the details 
of the parameters. One example is shown in Figure (7), where we compare the perturbation produced 
by mini-halos with that produced by the streams. The comparison is made for two values of ;n m h (red) 
and two values of A &t (blue). The intersection between blue and red lines shows the combination where 
one perturber type becomes dominant over to the other. In addition, in the case of having a smooth 
halo component, the fraction of mini-halos and the fraction of streams are no longer complementary, 
but the total perturbation effect can still be read from this plot. 

4 A possible relation between A st and 

The actual constraints for the streams like structures could be interpreted in terms of their original 
mini-halo mass if we had information about the disrupted time (t) and how the disruption was. These 
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Figure 7: Binary perturbations due to the presence of mini-halos and streams. For selected values 
of the stream linear mass density (A sU blue), an mini-halo mass (m m h, red), as a function of their 
corresponding sub-structure fraction. The intersection between blue and red lines shows the fraction 
where one perturber type becomes dominant over to the other. 



two aspects are not easy to determine. We propose a very simple approach to establish this connec- 
tion. Estimate the mass per unit length of a resultant stream, A st , by computing first its length L with 
the relation L ~ cr mh t [30], which assumes that the particles moves with the velocity dispersion of 
the progenitor mini-halo. For simplicity, we will consider the mini-halos had an initial NFW density 
profile, and so the maximum velocity dispersion can be inferred from the virial mass [59] through 

o* * 0.472V 2 ax ( 4.l) 
2 GM m 0.432 

max ry c 

2r -v Log[l + c]-- 

1 + C 

where r s , and r v ; r are the scale and virial radius, respectively. Use the following definitions: 
r v j r = (M v i r /(47r/3p cr Q m 5 t h)) 1 ^ 3 with p cr = 270.5Mq/Kpc 3 , the critical density of universe and 
<5 t h = 340 the predicted over-density of a collapsed object according to the top-hat collapse model 
for the A - CDM model with Q. m = 0.3); and, the concentration for mini-halos c = r v i T /r s = 12 * 
(M v i r //j _1 10 12 ) -012 [14]. Given the above relations, the resultant velocity dispersion for the halo 
mass in the range 10~ 6 < ra^M©) < 10° results to be 10 -4 < (j^ (km/s) < 10~ 2 , so after an 
evolution of a Hubble time, the length of these structures would be 7 < L (pc) < 370. Then, using a 
simple linear fit the relation between the mini-halo mass, and the linear mass density could be roughly 
set to 

A st * 10" 8 m^ 7 h . (4.2) 

For instance, an initially 1O~ 6 M0 mini-halo mass would be associated with a stream of linear density 
10 _13 Mq/AU, which according to Figure (3a) would be perfectly allowed by the dynamics in the 
Solar System. For comparison we have scaled our results in Figures (4), and (6a), using the relation 
in (4.2). However, in order to have a more accurate connection between A st and an statistical 
study for the interaction between mini-halos and perturbers should be performed, but this is out of 
the scope of this work. 
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5 Discusion and conclusions 



We have estimated the dynamical effect that sub-structure in the galactic halo may trigger on the evo- 
lution of the planet orbits in the Solar System; the Sun-Neptune and Earth-Moon orbits specifically. 
The effect has been compared to the current precision in the determination of their semi-major axis 
(known with millimeter accuracy for the Earth-Moon thanks to the APOLLO-LLR intiative [31]), 
finding that a local DM distribution completely populated by stream like structures with linear den- 
sity of A st > 10~ 10 Mq/AU in the solar neighborhood, can cause a secular evolution to the Earth-Moon 
system that would be above current measured limits (though Neptune orbit allows a higher value). We 
have also shown the allowed, and not allowed, parameter regions, in a "linear density-sub-structure 
fraction" space, for dark streams based on the current precision of the Earth-Moon distance measure- 
ments. 

Here it is worth to emphasize the importance of the local DM density value in our study. 
The typical approach to constrain its value is based on dynamical studies, however, many astro- 
physical uncertainties are associated. According to different authors the actual limits varies from 
« 10~ 18 M o /AU 3 ((0.2 - 0.8) GeVcm" 3 ) [43-47]. It is important, as well, to say that all these deter- 
minations only give information about the mean value, at the scale of 100's to 1000's of parsecs, but 
says nothing about the way the DM is distributed at smaller scales. Indeed, the density at the Solar 
System scales could show important over(under)-densities regardless the mean value at larger scales. 
Our study seems to indicate that if the Solar System DM were mainly distributed in mini-halos, the 
density could be an order of magnitude greater than the standard value without being in contradic- 
tion with the Solar System dynamics, but not if it is mainly distributed along streams. Also, note 
that as we have better kinematic constraints for stars in the solar neighborhood, the detection of dark 
sub-structure could be a reality; the GAIA astrometric mission can open that possibility [20, 21]. 

We have also applied our formalism to determine if the presence of Dark Substructure in dSphs 
Galaxies could be in conflict with the presence of wide open binaries inside dSphs. The case of 
multiple encounters, the diffusive regime, and the interaction with streams has not been considered 
elsewhere. We have demonstrated that catastrophic encounters are not, necessarily, the driven mech- 
anism of destruction of wide binaries. The effect of multiple encounters with small perturbers can 
cause secular evolution as well. Our results indicate that to have better constraints on the lifetime of 
the binaries could help to impose constraints to the properties of the dark sub-structure. But, as there 
are many parameter combinations, (m m h, A st , / s , t ev ), for which it is possible to cause secular evolu- 
tion of the binary system, or even dissolve it; we consider that the existence of this kind of binaries 
is not a sharp discriminant of the dark matter hypothesis. However, to get a good characterization 
of their distribution would be helpful to constrain the properties of the dark sub-structure, specially 
when combined with other astrophysical restrictions. 

Also, by relating the linear density of streams with a mini-halo progenitor mass (a simple ap- 
proach) we conclude that the typical cut-off in the mas power spectrum, i.e. the mass of the smallest 
structures ever formed in the Universe m m ; n ~ 10~ 4 - 10~ 12 Mq [9, 10], will correspond to streams 
with linear densities of A it ~ 10 -11 - 10~ 17 Mq/AU. Such values are compatible with the dynamics of 
the solar system (a similar conclusion was obtained when we consider the mini-halos as a point mass 
perturbers), and would be compatible also with the existence of wide binaries in dSph's, if they were 
to be detected. On the other hand, for streams with A « 10" 10 -10" 8 M©/AU (m^n ~ 10~ 3 - 10°M o ), 
the fractional energy imprinted to the binary stars would cause their disruption, and the corresponding 
effect to Earth-Moon system would be just below (or above) the observational restriction. That is, if 
wide binaries in dSph's with separations of a ~ O.lpc were not observed, it could be due to the pres- 
ence of streams with A st > 10~ 10 Mo/AU (m m i n ~ 10~ 3 ), but the presence of such kind of structures 
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in the solar neighborhood would (likely) be incompatible with the Earth-Moon dynamics (although 
it depends in the sub-structure fraction). This is an example of how the combined restrictions, from 
different astrophysical systems, could be revealing the sub-structure mass function cut-off scale. 

Finally, we want to make some comments regarding the direct, and indirect DM detection meth- 
ods. First, it is well known that the presence of dark mini-halos, or streams, in the solar neighborhood, 
could affect the event rate and energy spectrum, and consequently the exclusion limits set by current 
experiments. From our study, we can say that the dynamics of the Solar System is not significantly 
affected by the presence of both kind of sub-structures, mini-halos and streams (not even in the ide- 
alized case when the sub-structure fraction is equal to one). This rise up the possibility of being in a 
over-dense region and suggest that DM halo models that includes the presence of DM sub-structure 
should be used when interpreting experimental results. However, there is also the possibility to be 
in a under-dense region, not distinguishable from the standard case by dynamical studies, but with 
important consequences for the direct detection methods. 

On the other hand, the presence of dark sub-structure in galactic halos, could also affect the 
interpretation of the results for DM indirect detection experiments, because in some cases it is sub- 
jected to the use of a boost factor mechanism, usually associated to the Sommerfeld enhancement, 
that is velocity dependent. 3 As mini-halos and streams are, in general, kinematically colder than a 
soft galactic halo, then, the limits to the sub-structure properties set in this work could be useful also 
to interpret the results of DM indirect detection experiments . 
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A Corrections for Spherical models. 
A.l Correction for extended perturbers. 

In this section we present the behavior of the U function, in equation (2.1), that accounts for the 
extension and structure of the dark mini-halos. For this we have assumed the mini-halos have a NFW 
mass profile, normalized to the virial mass, of the form: 



where c = r va /r s is the concentration; r s and r v i r are the scale and virial radius respectively. Fol- 
lowing [35], and given that the mass profile is normalized to the virial mass, we will define the U 
function as: 
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Figure 8: Structure function, U, for halos with a NFW density profile, as a function of their mass. 
For comparison we added a line for the point mass case (£/nfw = 1), to see the scale at which this 
becomes valid. 

After integration, the U function takes the form 



and it is ploted in Figure, 8, for different values of the mini-halo mass. The U function crosses the 
point mass approximation at r/p = c for the different masses. For comparison we added a line for 
the point mass case (£/nfw = 1), to see the scale at which this becomes valid. We have seen this 
correction is less important for small mini-halos. 

A.2 Corrections for non-distant tide encounters. 

In this section we deal with encounters where the distant tide approach is not applicable. In particular, 
this happens when we study the interaction of wide binaries with point mass (or spherical) dark 
sub-structure present in the MW satellites. We will use the this system to state the problem, and 
the solution, but this can be applied to other systems for which the distant approximation does not 
applies. 

The distant tide approximation is valid when the size of the target system is much smaller than 
the impact parameter. Te determine whether this condition is satisfied, or not, we consider the char- 
acteristic of the wide binaries and DM sub-suucture we used in Sec.3.2 (p dm = 2.2Kr 17 M /AU 3 , 
a=0.lpc, Mj, = IMq, etc.). We can compare the binary separation with the mean impact parameter 
for the encounters, < p >= T(l/3) (367rpdm/"imh/ 3 ^ (defined from the nearest neighbor distribu- 
tion), for different values of the perturber mass. Table 9a shows this comparison: values lower or 
close to one implies that distant tide approximation is not applicable. For the less massive perturbers, 
most of the encounters could not be treated in the distant tide approximation. Since the sampling 
of impact parameters is random, we do not know, a priori, which of the encounters could be treated 
in the distant approximation and which not. In the practice we have to make this comparison for 
each encounter, then we decide if it is appropriate to use the distant tide approximation, or to directly 




for p/r s < c 



(-£ + log[l+c]) 



, 1 for p/r s > c, 



(A.3a) 
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calculate the velocity kick to each component of the binary, and calculate the input energy as we 
describe in the following paragraphs. 




(a) (b) 



Figure 9: (a) Shows the factor p me!m /a for different values of the mass perturber. Values close to one 
means that the distant approximation, in the calculation of the dynamical effect of encounters between 
a target and a population of mimi-halos, is not applicable, (b) Diagram that shows the definitions used 
to calculate the correction for non-distant encounters with mini-halos. 



For each component of the binary the kick in velocity due to an encounter, in the non-distant 
tide approach, is given by 

2GM„ p. 

Av; = P -U{pd^ (A.4) 

V Pi 

where U(p) is a function that depends on the perturber density profile, and tends to one when the 
impact parameter (p) is larger compared with the size of the perturber (this correction is described in 
detail in [35], and in appendix A.l). The subscript i refers to the i-th component of the binary. The 
differential energy per unit mass transferred to the binary system is 

AE = 0.5Av^ iff + v • Av (A.5) 

where Avdiff = |Av2 - Avi |, is the vector difference of the individual contributions of both components 
of the binary system. 

In our approach the magnitude of the impact parameter is sampled from a nearest neighbor 
distribution, so we have to assume respect to what point point of the binary it is referred, by simplicity 
we will assume it is respect to the center of mass, so that the individual impact parameters to the 
binary components are given by: 



p\ = p 2 + (a/2) 2 - pa cost] 

p\ = p 2 + (a/2) 2 + pa cost] (A.6) 

where 77 is the relative direction between the imaginary line that connects the binary, and the one that 
draws the impact parameter direction (we will sample this angle from a uniform distribution between 
and 7r, since it is equally probably to have any orientation of the binary). Given this geometry we 
rewrite eq. (A.5) in terms of the impact parameters and the binary separation. 
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AE = 0.5 



(2GM„\ 2 , 

I— —I <5vf ot where 



$ v tot = & v \ + ^ v 2 ~~ 2(5Vi • 5V2 

^l)\ 2 , /^2)\ 2 2 ^/(Pl)^(P2)(^+^-^ 2 ) 
Pi / \ PZ I p\p\ 

This last expression is used whenever the condition of distant tide encounters, in our analysis, 
is not satisfied. Note we have neglected the linear velocity term in eq (A.5), beacuse we expect it 
approaches to zero when AE is summed over all the encounters as a result of symmetry (See [35]). 

B Stream models 

B.l Derivation of tidal perturbation. 

In this appendix we derive analytical expressions for the force and tidal fields produced by the 
stream models we have considered. All of them present axial symmetry. Our simplest model is a 
1-dimensional mass distribution with linear mass density A st (S-1D). Then we consider a stream with 
finite cross section of radius R and constant cross-sectional density S (S-CD). Finally, we consider 
streams whose cross-sectional densities vary as power laws with central cores (S-Core), or cusps 
(S-Cusp): 

Z(R;R ,a,y) = E (R/R o y a [1 + (/?//? )f~ r , (B.l) 

where S is a scaling density parameter, R Q a characteristic radius and a, y are dimensionless expo- 
nents that satisfy: 

< a < 2, y > 2. (B.2) 

The restriction a > is to avoid an inverted density profile at the center, while a < 2 avoids an 
infinite mass singularity there (0 < a < 2 produces a finite mass singularity). The restriction y > 2 is 
to have a finite mass per unit length. When a = we have the flat core cases, whereas for a > we 
have the cusp cases. 

Force field 

The force field can be obtained using the integral Poisson equation: 



V -FdV = -AnGM. (B.3) 



Taking advantage of the axial symmetry of the problem, the resulting forces are given by: 

F s ~ lD (R)= (B.4) 

F s-c„„ (S) . _ Fc „--r,(i,,) + (i + ,F-i 

77(1 +7]Y 

F s-Cus P(R)= _ Fo 1 V_ (a=1)? (B _ 7) 
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R/R 

(a) Magnitude of the force produced by our stream 
models as a function of distance to the stream cen- 
tral axis, R. It coincides with the tidal force in the 9 
direction, equation (B.10). 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 



(b) Magnitude of the tidal force, in the R direction, 
produced by our stream models as a function of dis- 
tance to the stream central axis, R. 



Figure 10: The curves are for the 1-D (black), constant density (blue), core (green) and cusp (red) 
models. In each case, we set A st , F$ and £o to the unity. For the core and cusp models y = 3 is 
assumed. 



where we have introduced the dimensionless length r\ - R/R () and F Q = 2GA st /R , with A st equal to 
the total mass per unit length (integrated across the stream cross section). All forces are orthogonal 
to the stream and directed towards it. For the cusp case there is no general analytical expression, so 
we show the special a = 1 case, which is the only cusp model we will consider from now on. 

Figure (10a) shows these forces, for streams with A st = 1, or Fq = 1, accordingly; and y = 3 for 
the core and cusp models is assumed. We see that the S-1D and S-CD models coincide for/? > R a , but 
whereas the former diverges toward the center, the latter converges linearly to zero. The core model 
goes also to zero at the center, while the cusp model does not and so the force is discontinuous at the 
stream axis. For large R all models coincide, as they should, since they become indistinguishable. 

Tidal force field 

The components of the tidal force produced by a force field F on a target of size dr at position r, 
in the linear approximation, (we use general coordinates x a and Einstein's sumation convention), is 
given by : 

dF 

F, Xa = -^dyP = T afS dxP (B.8) 

The tensor t, a symmetric and second rank tensor, is the Jacobian of the gravitational force field F, 
and the negative of the Hessian matrix, H, of the scalar potential function O. 4 The tensor r can 
be written in any curvilinear coordinates by means of the covariant derivative, [See appendix A of 
reference 62, for a detailed explanation], i.e. 

Tafi = -H aP «D) = F Qfi - T y ap F a , (B.9) 

where, F = -V®. The Christophel symbols, T, and the force F, must be written in the selected 
coordinate system. For our specific case of streams with axial symmetry, it is convenient to use 
cylindrical coordinates, defined as in Figure (11a); then, in this particular case the tidal force takes 
the form: 

4 The vector r spans the range across where the differential effect of the tide is to be calculated. 
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F t (r,dr) = T rr dre r + TggrdOeg 
dF{f) 



dr 



dre r + F(r)d9e e . (B.10) 



Using the expressions for the force of each of our stream models, eqs. (B.4) to (B.7) respec- 
tively, it is straightforward to obtain the expressions for the corresponding tidal forces, given by eq. 
(B.10). As an example, the tidal force in the 1-D stream case is, 

rS-iD 2G/l s t A 2G/l s t m 
Ff lu = — ^— dre r —d6t e , (B.ll) 

Figure (10b) shows radial component of the tidal forces, for streams with A st - 1, or Fq - 1, 
accordingly; and y - 3 for the core and cusp models is assumed. For the S-1D and S-CD models is 
the same stretching force outside the stream; once inside the finite thickness stream, the tidal effect 
turns into a constant compression. The Core stream presents a simmilar behaivor as the S-CD, i.e, in 
some range gives a compresive force, and in other a streaching one. The Cusp stream gives always a 
compresive force. The tangential component of the tidal force coincides with the force, and the effect 
is always compressive for all the models. 

The fractional energy change 

Now, to quantify the effect of an encounter with a stream we will focus on changes to the internal 
energy of the target, which is mostly kinetic. One needs to compute the integral of the tidal force over 
time to get the kick in the velocity during to the encounter with a stream, 



<Sv = J F,(R,<5r)dt, (B.12) 

where R is the cylindrical position vector of the target center with respect to the closest stream point, 
and Sr, is the position vector of an element of the target, with respect to the target center. Both, R and 
6r are functions of time. The integral in eq.(B.12) should be made over the time that the encounter 
takes place, in general this could be accomplished by means of numerical simulations. However, this 
problem can be treated analytically in two opposite regimes: If the encounter time is a lot shorter than 
the internal dynamical time of the target, then we can take 6r as fixed and just take into account the 
variation of R(/), this is the so called, impulsive regime. The opposite regime is when the encounter 
time is a lot longer than the internal dynamical time of the target. In this case we can average the 
effect of the tidal force over an entire dynamical timescale of the target and then use the average force 
as the integrand in eq. (B.12), this is the adiabatic regime. 

The systems of interest to this work are suitable to be worked in the impulsive regime, so we 
will estimate the effect of a single encounter using this approximation. The kick in velocity will 
be obtained by eq. (B.12) once we substitute the corresponding tidal force, (B.10), for the different 
models. But before doing this we must consider the geometry of the problem. There are two issues 
here: first, the geometry of the encounter stream-target will define the form of R(f); and second, the 
orientation of the target with respect to the stream will define the components of dr. 

The geometry of the encounter stream-target Let us introduce here the concept of the impact 
parameter, p, for encounters with streams. It will be defined as the closest distance between the 
stream and the target, with the particularity that it will always be orthogonal to the to the z axis of the 
cylinder. Now, lets see the relation of the impact parameter with the vector position R of the target 
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(b) 



Figure 11: (a)Coordinate frame used for the tidal tensor computation, (b) Relative direction of 
movement definition. 

with respect to the stream. Lets suppose we define the impact parameter at the time t - ti, it should 
lie in the plane normal to z. At this time the impact parameter, p, and the vector position R coincide. 
Now as the time evolve, the target will move, in some direction with respect to the stream, so that 
the vector position will not coincide any more with the impact parameter. The difference between 
them will be just the distance the target, moving at velocity vo relative to the stream, has traveled. 
This distance will depend the relative direction of movement between the stream and target, being 6 
the angle that characterize this relative direction, as it is shown in Figure (lib). So, the equation that 
gives the temporal variation of the target-stream magnitude vector position, R, is: 



One can verify that 8-0 correspond to the case in which the stream and target are moving parallel, 
at constant impact parameter. 

The geometry of the target orientation Now it is time to turn our attention to the orientation of the 
target with respect to the stream. What we need to do is to express the vector position of an element 
of the target, Sr, in the reference frame defined by the stream (Figure 11a), that is to express 5r with 
respect to the triad of unit vectors (e^, eg, e z ). 

One can describe this transformation by means of a three angle rotation: 

1. Rotation of the target plane by an angle a with respect to the eg axis (Figure 12a). 

2. Rotation of the target plane by an angle /? with respect to the original e z axis (Figure 12b). 

3. Rotation on the target plane by an angle \]j with respect to the normal e n of the target plane 
(Figure 12c). 

The first 2 angles defines the spatial orientation of the target plane with respect to the stream 
plane, the third angle is just the orbital phase of the target, see Figure 12 for clarification in the 
definition of the angles. It is clear that -7r/2 < a < n/2, < ft < n, and < ifr < 2n. The rotation 
matrix for the individual rotations are: 




(B.13) 
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( cos OS) sin OS) 0^ 
-sin 05) cos OS) 
1 



1 cos (a) - sin (a) > 
1 
sin (a) cos (a) 



cos(i/f) sin(i/r) 0' 
- sin (if/) cos (if/) 
1 



(B.14) 



To get the proper final rotation, we should apply them in the order: first the yS-rotation, followed 
by the c-rotation and finally the (/^-rotation. The inversion in the order of the first two is so that the 
pivot for the a-rotation is with respect to the rotated eg axis. If we do it in the inverse order, then 
the pivot for the second rotation OS) is the original e z axis, and not its transformed version, which 
complicates the algebra. The expression of vector 6r using in reference frame attached to stream, 
then, is: 



6r = 8r {[cos (a) cos OS) cos (if/) - sin OS) sin (if/)] e« 

+ [cos (a) cos OS) sin (if/) + cos (/?) sin (if/)] eg - [cos (if/) sin (a)] e,} 



(B.15) 




Figure 12: We show the definitions of the rotation angles a, ft and if/ that allow us to write the vector 
position of an element of the target, 8r (that lies in the target plane -green-), in the stream reference 
frame coordinates (brown plane). 



Finally, the perturbation energy (per unit mass) transfered to the internal dynamics of an ex- 
tended target system, by one single stream encounter, is simply 

AE = A(l/2v 2 ) = (v • Av) + l/2(Av) 2 ~ l/2(Av) 2 (B.16) 

where in the last step, we have neglected the linear term between the velocity of the target 
planet, v, and the tidal velocity kick, Av (random orientation). We now compare this perturbation 
energy with the binding energy of the target: E\, = -G M c /(2a), and substitute equations (B.15) and 
(B.10) in eq.(B.12). The integration was done considering 6r constant, i.e. dr,a,fi and if/ constants 
during the encounter. The results are the following, for the different stream models. 



(Aff) st 
\E b \ 



4Gn z A 



2)2 



M c v 2 srn(0)V 



< T(a,ifr) (ID) 

T(a,if/)B 2 (R fp) (CD) 

C(R Q /p,a,p,if/) (Core) 

{ B(Ro/p,a,f3,if/) (Cusp). 



(B.17a) 
(B.17b) 
(B.17c) 
(B.17d) 



where the geometric functions, T, B, C, and D are given by: 
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P/Ro P 

(a) Perturbation energy produced by the streams, as a (b) Impact parameter probability distribution functions 
function of the impact parameter, p. The models are for dark streams calculated by the sampling procedure 
the same as in Figure 10a described in the text. 

Figure 13 
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(B.18) 
(B.19) 

(B.20) 



(B.21) 



All the functions AE/Ej, for the different models, converge to that of the S-1D model for large 
impact parameters, see Figure 13a. 

B.2 Impact parameters for streams. 

In order to generate the distribution function of impact parameters for the streams, we follow the next 
procedure. We start by tossing a given number of random lines within a region of the 3D-space. Then 
we also tose a collection of uniformly distributed points, within the same region, and determine the 
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distance from each point to the closest line, that is the impact parameter. This results in a sampling of 
the underlying distribution of impact parameters that is completely independent of the introduction 
of a sample volume. We show the resultant probability distribution in Figure (13b). We still need to 
relate this realization to a specific distribution of matter with a certain spatial mass density, p, and 
linear massdensity for the lines, A, so that the impact parameter denoted by p* can have the correct 
dimensions according to the specific physical situation. The way to do this to introduce a sample 
volume to measure the mass density around each of the tests points. We did several experiments, 
keeping fixed the total number of lines in the realization, as well as the realization volume, but 
changing the size of the sample volume, R*, and we found that the number density of streams varies 
with it, but the mass density doe not. Thus, in our procedure, the mass density is a well defined 
physical quantity whose value is independent of the size of the region used to compute it. In the 
practice, we normalized the linear mass density of the streams to the unity so that the mass density 
we get from the realization is p*. By doing a dimensional analysis we get to the conclusion that 
the relation between the impact parameter, p* , and the physical impact parameter p, is given by: 
P = ( y^stPVfnnsp) P*> where A st and p are, given in physical units, and p* is the numerical mass 
density obtained in a given realization. 

We will use the impact parameter distribution function shown in Figure 13b, that corresponds to a 
numerical mass density of p* = 1205.08; the A st , f^, and p parameters will be fixed according to the 
specific situation under study. 
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